GNSS post positioning with elongated dither sequence

ABSTRACT

A computer apparatus for post positioning with a selected precision. The apparatus includes a GNSS post processor to post process reference GNSS carrier phases from a reference system and rover GNSS carrier phases from a rover receiver to compute a secure position for the rover receiver not available to a user. The apparatus includes a random process generator to generate a sequence of offset vectors to dither the secure position according to a computed dither level to provide the selected precision for a user-available position for the rover receiver.

CROSS REFERENCE TO RELATED APPLICATIONS

This application is a continuation in part of pending application Ser. No. 11/799,022 by Bird filed Apr. 30, 2007 which is a continuation of application Ser. No. 11/138,223 by Bird filed May 26, 2005 now patent U.S. Pat. No. 7,227,496 issued Jun. 5, 2007, all assigned to the same assignee.

BACKGROUND

1. Field of the Invention

The invention relates generally to positioning and more particularly to positioning with selected accuracy having high integrity.

2. Description of the Background Art

The Global Positioning System (GPS) is operated by the United States government for providing free GPS positioning signals to all users around the world. Stand alone GPS receivers can use a coarse/acquisition (C/A) code in these signals for computing unaided positions having typical accuracies of about five to twenty meters. These accuracies are sufficient for some applications including most navigation applications. However, there are positioning applications, such as survey, mapping, machine control and agriculture, where greater accuracy or integrity is needed.

Some of these needs are met by differential GPS systems that provide GPS code phase corrections. A GPS receiver that is constructed for differential GPS operation can use the code phase corrections for computing positions having typical accuracies of a few tens of centimeters to a few meters. These accuracies are sufficient for many positioning applications. However, a user cannot be altogether confident in the accuracies of stand alone or differential GPS positions because the integrity of the positions is affected by multipath. Multipath reflections of the GPS signals can cause occasional large errors of tens to hundreds of meters or even more depending on the extra distances that are traveled by reflected signals.

Fixed ambiguity real time kinematic (RTK) systems provide highly accurate GPS carrier phase measurements in order to provide greater accuracy and at the same time avoid most of the effects of multipath. A rover GPS receiver that is constructed for RTK operation can use the carrier phase measurements for determining relative positions having typical accuracies of about a centimeter to a few tens of centimeters. The term “fixed ambiguity” refers to the fact that an integer number of cycles of carrier phase is resolved (fixed) for the RTK carrier phase measurements between the reference phase and the phase measured by the rover. The resolution of the carrier cycle integer traps multipath signal errors that are greater than a portion of the wavelength of the carrier of the GPS signal, resulting in a high confidence and integrity for the RTK-based positions.

Existing GPS RTK systems provide fixed RTK carrier phase measurements to the users for a cost that is largely driven by the fixed infrastructure costs for providing the system divided by the number of users. However, some users require the integrity of fixed RTK-based positioning but do not require the full accuracy that it provides. Unfortunately, there is no existing technique for spreading the infrastructure costs across more users by providing high integrity positions with accuracies that are lower than the full accuracy of the system.

SUMMARY

The present disclosure describes ways of providing high integrity positioning with controlled accuracies for a rover station either by providing synthetic reference phases for a GPS reference system or by dithering a secure rover position.

Briefly, several systems are disclosed using measurements of or including one or more real time kinematic (RTK) reference stations for receiving GPS signals at one or more actual reference positions and for measuring reference phases. When three or more reference stations are used, virtual reference phases may be determined for a virtual reference position. A synthetic offset vector is generated in a reference station, a server in the reference system, an RTK rover station, or a synthetic phase processor interacting between the reference stations and the rover station. Reference phase measurements are used with the synthetic offset vector for inferring synthetic reference phases for a synthetic position where the synthetic position is not equal to any of the actual or virtual reference positions. The rover station uses the actual or virtual reference position with the synthetic reference phases in place of the actual or virtual reference phases for computing a rover position with respect to the actual or virtual reference position having an added positional error that is proportional to the synthetic offset vector.

In another approach a secure RTK rover station uses a synthetic offset vector directly for dithering a secure rover position determined from the actual or virtual reference phase. The synthetic offset vector may be generated in a reference station, a server in the reference system, the rover station, or a processor acting between the reference system and the rover station. The positions determined by the rover station have the integrity of the RTK system with accuracy controlled by the synthetic offset vector.

One embodiment is a secure rover station having a controlled accuracy for a geographical position, comprising: a rover global navigation satellite system (GNSS) receiver for determining a secure position not available to a user of the rover station; and a position dither processor for dithering the secure position with a selected non-zero synthetic offset vector for issuing a rover position available to the user having an added position error proportional to the synthetic offset vector.

Another embodiment is a method for controlling accuracy of a geographical position, comprising: receiving a global navigation satellite system (GNSS) signal; using the GNSS signal for determining a secure position not available to a user of the rover station; and dithering the secure position with a selected non-zero synthetic offset vector for providing a rover position having an added position error proportional to the synthetic offset vector to the user.

Another embodiment is a tangible medium containing a set of instructions for causing a processor to carry out the following steps for controlling accuracy of a geographical position, comprising: receiving a global navigation satellite system (GNSS) signal; using the GNSS signal for determining a secure position not available to a user of the rover station; and dithering the secure position with a selected non-zero synthetic offset vector for providing a rover position available to the user having an added position error proportional to the synthetic offset vector.

Another embodiment is a computer apparatus for post positioning with a selected precision, comprising: a GNSS post processor to post process reference GNSS carrier phases from a reference system and rover GNSS carrier phases from a rover receiver to compute a secure position for the rover receiver not available to a user; and a random process generator to generate a sequence of offset vectors to dither the secure position according to a computed dither level to provide the selected precision for a user-available position for the rover receiver.

Another embodiment is a method for post positioning with a selected precision, comprising: post processing reference GNSS carrier phases from a reference system and rover GNSS carrier phases from a rover receiver for computing a secure position for the rover receiver, the secure position not available to a user; and generating a sequence of offset vectors to dither the secure position according to a computed dither level for providing the selected precision for a user-available position for the rover receiver.

Another embodiment is a computer-readable medium having computer-executable instructions stored or carried thereby that when executed by a processor perform a method comprising steps of: post processing reference GNSS carrier phases from a reference system and rover GNSS carrier phases from a rover receiver for computing a secure position for the rover receiver, the secure position not available to a user; and generating a sequence of offset vectors to dither the secure position according to a computed dither level for providing the selected precision for a user-available position for the rover receiver.

These and other embodiments and benefits of the present invention will no doubt become obvious to those of ordinary skill in the art after having read the following best mode for carrying out the invention and viewing the various drawings.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 is a diagram of a single reference system of the prior art for providing reference phases to a rover station;

FIG. 1A is a diagram of a reference network positioning system of the prior art for providing reference phases to a rover station;

FIG. 2 is a diagram of a single reference positioning system for providing synthetic reference phases to a rover station for adding a position error to a rover position;

FIG. 2A is a block diagram of a reference station for the system of FIG. 2;

FIG. 3 is a diagram of a reference network positioning system for providing synthetic reference phases to a rover station for adding a position error to a rover position;

FIGS. 3A, 3B and 3C are block diagrams of first, second and third embodiments of a server for the positioning system of FIG. 3;

FIG. 4 is a diagram showing a secure rover station for computing synthetic reference phases for adding a position error to a rover position operating in a single reference positioning system;

FIG. 4A is a block diagram of the rover station of FIG. 4;

FIG. 5 is a diagram showing a secure rover station for computing synthetic reference phases for adding a position error to a rover position in a reference network positioning system;

FIGS. 5A and 5B are block diagrams of first and second embodiments of the rover station of FIG. 5;

FIGS. 6A and 6B are block diagrams of first and second embodiments of random reference generators for providing synthetic reference phases;

FIG. 7 is a diagram for a single reference positioning system where a secure rover station dithers a secure position for providing an unsecure position having an added positional error;

FIGS. 7A and 7B are block diagrams of first and second embodiments of the positioning system of FIG. 7;

FIG. 8 is a diagram for a reference network positioning system where a secure rover station dithers a secure position for providing an unsecure position having an added position error;

FIGS. 8A, 8B, 8C and 8D are block diagrams of first, second, third and fourth embodiments of the positioning system of FIG. 8;

FIG. 9 is a block diagram of a random position dither processor for the systems of FIGS. 7 and 8;

FIG. 10 is a flow chart of a method for providing synthetic reference phases from a single reference positioning system to a rover station;

FIG. 11 is a flow chart of a method for providing synthetic reference phases from a reference network positioning system to a rover station;

FIG. 12 is a flow chart of a method for computing synthetic reference phases in a rover station operating in a single reference positioning system;

FIG. 13 is a flow chart of a method for computing synthetic reference phases in a rover station operating in a reference network positioning system;

FIG. 14 is a flow chart of the method for dithering a secure rover position for providing an added error to a rover position in a single reference positioning system;

FIG. 15 is a flow chart of the method for dithering a secure rover position for providing an added error to a rover position in a reference network positioning system;

FIG. 16 is a block diagram of a GNSS positioning system having a reference system, a rover receiver, and a post processing computer apparatus;

FIG. 17 is a position diagram showing secure and user-available rover positions for the rover receiver of FIG. 16;

FIG. 18 is a block diagram of the computer apparatus of FIG. 16;

FIG. 19 is a block diagram of a random process generator of the computer apparatus of FIG. 18;

FIG. 20 is a block diagram of a vector offset generator and a position dither processor of the computer apparatus of FIG. 18;

FIG. 21 is a block diagram of a DOP scaler for the vector offset generator of FIG. 20;

FIG. 22 is a flow chart of a method for providing a rover position having a selected precision that is available to a user;

FIG. 23 is a flow chart of a method for generating an elongated sequence of offset vectors for dithering a secure rover position.

FIG. 24 is a flow chart of a method for dithering a secure position for providing a user-available position; and

FIGS. 25A and 25B are flow charts of a method for selecting a user-available precision based on a DOP of a secure position.

DETAILED DESCRIPTION

The details of several embodiments for carrying out the idea of the invention will now be described. It should be understood that the description of these details is not intended to limit the invention to these details. On the contrary these details are merely intended to describe the best mode for carrying out the idea of the invention. Numerous alternatives, modifications and equivalents of the embodiments described herein will be apparent to someone skilled in the art as within the scope of the idea of this invention. An embodiment of the invention is described for the global positioning system (GPS). However, it will be apparent to those in the art that the invention may be carried out with a generic global navigation satellite system (GNSS) including the global positioning system (GPS), the global orbiting navigation system (GLONASS), the Galileo system or a combination of these systems. It should also be noted that pseudolites may be used in place of satellites for broadcasting GNSS positioning signals.

FIG. 1 is a diagram showing a conventional real time kinematic (RTK) global positioning system (GPS)-based system of the prior art. The reference station 12 includes a reference GPS receiver for receiving GPS signals 14, illustrated as 14A, 14B and 14C, from GPS satellites 16, illustrated as 16A, 16B and 16C. The reference station 12 measures the carrier phases of the GPS signals 14 and sends a radio signal 17 having reference data for the measured phases and reference geographical position to one or more rover stations shown as rover 18. The rover station 18 includes an RTK GPS receiver for measuring the carrier phases for the same GPS signals 14.

The difference between the reference and rover phase measurements yields estimates of perpendicular distance vectors represented by a vector d between the rover station 18 and the GPS satellite 16A. Measurements from several GPS satellites 16 yield estimates of several perpendicular distance vectors and ultimately the position of the rover station 18 with respect to the reference station 12. The vector d may be understood as a dot product of the vector between the GPS satellite 16 and the rover station 18 and the vector between the reference station 12 and the rover station 18. An exemplary RTK GPS system is described in U.S. Pat. No. 5,519,620, entitled “centimeter accurate global positioning system receiver for on-the-fly real-time-kinematic measurement and control” by Nicholas C. Talbot et al., incorporated herein by reference.

FIG. 1A is a diagram showing a virtual reference system (VRS) RTK GPS-based system of the prior art. Reference network stations 12A, 12B through 12N include reference GPS receivers for measuring carrier phases of the GPS signals 14 received from the GPS satellites 16. The reference stations 12A-N send signals 22 having reference data for their measured phases and reference geographical positions to a server 23. One of the reference stations, illustrated as 12A, is designated as a master reference station. The server 23 and the master reference station 12A may be located together. The server 23 communicates with one or more VRS RTK GPS rover stations with a radio signal 25. The VRS RTK GPS rover stations are shown as rover 24.

The server 23, or the server 23 together with the rover station 24, determine a virtual reference position 26 and a virtual vector 27 between the position of the master reference station 12A and the virtual reference position 26; and then uses the position and measured phases of the master reference station 12A, the positions and measured phases of the auxiliary reference stations 12B-N and the virtual vector 27 (or virtual reference position 26) to calculate virtual reference phases for the virtual reference position 26 according to a virtual reference system (VRS) parametric model. The rover station 24 includes an RTK GPS receiver for measuring the phases for the same GPS signals 14. The difference between the virtual reference and rover phase measurements yields estimates of perpendicular distance vectors to the GPS satellites 16 analogous to the vector d, described above, and ultimately the position of the rover station 24 with respect to the virtual reference position 26.

The use of a network of reference stations instead of a single reference allows modeling of the systematic ionosphere and troposphere parametric errors in a region and thus provides the possibility of error reduction. Networks exist using public domain RTCM and CMR standards for bi-directional communication reference data to the rovers. Detailed information on the modeling of the errors is available in “Virtual Reference Station Systems” by Landau et al., published by the Journal of Global Positioning Systems for 2002, Vol. 1, No. 2 pages 137-143.

FIG. 2 is a diagram showing a real time kinematic (RTK) GPS-based positioning system referred to with a reference number 30. The positioning system 30 includes at least one reference station 31 for receiving the GPS signals 14 from the GPS satellites 16. The reference station 31 has a reference position that is established by a survey or some other means. The system 30 receives or generates a synthetic offset vector 32 for controlling the positioning accuracy that is provided by the system 30 for one or more RTK GPS rover stations 118. The reference position and the synthetic offset vector 32 define a synthetic position 33 where the synthetic position 33 is separated from the reference position by the synthetic offset vector 32. The length and direction of the synthetic offset vector 32 are arbitrary. However, it is normally a few meters or less.

FIG. 2A is a block diagram showing the reference station 31 and the rover station 118. The reference station 31 includes a reference GPS receiver 34, a synthetic phase processor 35 and a radio transceiver 36. The reference GPS receiver 34 and the synthetic phase processor 35 may be separate or combined into a single unit. The radio transceiver 36 may be a transmitter without a receiver if two-way transmission with the rover station 118 is not required. The reference GPS receiver 34 measures the carrier phases of the GPS signals 14. The processor 35 uses the synthetic offset vector 32 with the reference position for the station 31 and the three dimensional angles to the GPS satellites 16 for inferring the synthetic reference phases for the carrier phases that would be measured if the measurements were made at the synthetic position 33.

The radio transceiver 36 sends a radio signal 37 having synthesized reference data for the synthetic reference phases and the reference position to the rover station 118. Cellular or landline telephones may be used to provide or to augment the radio signal 37. The synthesized reference data still includes a correct geographical reference position of the reference station 31, as is conventional, but the phases for the GPS signals 14 are not the actual phases that are measured at the reference position but are instead the synthetic reference phases that are calculated from the measured reference phases, the actual reference position and the synthetic offset vector 32 (or synthetic position 33).

The rover station 118 includes a rover GPS receiver 119 and an anomaly detector 120. The rover GPS receiver 119 receives the GPS signals 14 and measures the carrier phases from the same GPS satellites 16 and computes the differences between the measured rover phases and the synthetic reference phases. Using the synthetic reference phases in place of the actual reference phases, it now arrives at estimated perpendicular distance vectors represented by the vector d* instead of the estimated perpendicular distance vectors represented by the vector d described above. The vector d* may be understood as a dot product of the vector between the GPS satellite 16A and the rover station 118 and the vector between the synthetic position 33 and the rover station 118. When the rover station 118 calculates its position with respect to the reference station 31, it arrives at a position 38 that has an added positional error 39 of equal length in the opposite direction as the synthetic offset vector 32. Using this technique the positioning system 30 is able to arbitrarily introduce the added error 39 into the position 38 that is calculated by the rover station 118.

The rover GPS receiver 119 determines double difference phase residuals from current and previous synthetic reference phases and current and previous measured rover phases and passes the phase residuals to the anomaly detector 120. The anomaly detector 120 detects a phase residual anomaly when the phase residual is greater than a phase threshold corresponding to a selected distance for an integrity limit 40 for RTK operation. The integrity limit 40 corresponds to an outer limit of a zone about the rover position 38. When an anomaly is detected, the anomaly detector 120 inhibits the rover GPS receiver 119 from providing the rover position 38 to the user of the rover station 118, or provides a notification to the user that the anomaly was detected and allows the user to decide whether or not to use the position 38. Alternatively, the anomaly detector 120 provides a solution for the rover position 38 where the synthetic reference phase and measured rover phase for the particular GPS signal 14 associated with the anomaly are not used. The effect of the system 30 is that the rover position 38 has the controlled added position error 39 without degrading the integrity limit 40 of the RTK positioning solution of the rover position 38.

FIG. 3 is a diagram showing a network embodiment of a real time kinematic (RTK) GPS-based positioning system referred to by a reference number 50. The positioning system 50 includes a network of reference stations, referred to as 51A, 51B through 51N, for receiving the GPS signals 14 from the GPS satellites 16. The reference network stations 51A-N have reference positions established by a survey or some other means. The system 50 receives or generates a synthetic offset vector 32 for controlling the positioning accuracy that is provided by the system 50 for one or more RTK GPS rover stations shown as a rover station 124A, 124B or 124C. The rover station 124A, 124B or 124C includes a rover GPS receiver 125A, 125B or 125C, respectively, and an anomaly detector 126A, 126B or 126C, respectively.

The positioning system 50 also includes a server 52A, 52B or 52C. The server 52A-C and the reference network stations 51A-N communicate with radio signals 54. One of the network stations, for example the station 51A, may be designed as a master and the other reference network stations 51B-N as auxiliaries. The master reference station 51A and the server 52A-C may be co-located and may or may not share processing power; or the master reference station 51A and the server 52A-C may be physically separated.

The reference network stations 51A-N measure carrier phases of the GPS signals 14 and then communicate their phase measurements to the server 52A-C. The server 52A-C communicates with the rover station 124A-C with a radio signal 56. In a conventional system, the server 23 uses the virtual vector 27 and the master and auxiliary reference positions and phases for determining the virtual reference phases for the virtual reference position 26. The sum of the virtual vector 27 and the synthetic offset vector 32 is a master synthetic vector 64. The position of the master reference station 51A and the master synthetic vector 64 define a synthetic position 133. The system 50 uses the synthetic offset vector 32 and the virtual vector 27 (or the master synthetic vector 64) and the master and auxiliary reference positions and phases for determining the synthetic reference phases for the synthetic position 133.

The rover station 124A-C expects reference phases as if the phases are measured at the virtual reference position 26, however, it receives the synthetic reference phases inferred for the synthetic position 133. The rover GPS receiver 125A-C receives the GPS signals 14 and measures the carrier phases from the same GPS satellites 16; and computes the differences between the measured rover phases and the synthetic reference phases. Using the synthetic reference phases in place of the virtual reference phases, it now arrives at estimated perpendicular distance vectors represented by the vector d* as described above. When the rover station 124A-C calculates its position with respect to the virtual position 26, it arrives at the position 38 relative to the virtual reference position 26 that has the added positional error 39 that is equal in length and in the opposite direction as the synthetic offset vector 32. Using this technique the positioning system 50 is able to arbitrarily introduce the added error 39 into the position 38 that is calculated by the rover station 124A-C.

The server 52A-C and the rover station 124A-C may use two-way communication to agree on the geographical position for the virtual reference position 26. For example, the virtual reference position 26 may be selected to be the best estimated position of the rover station 124A-C. It should be noted that the invention is not dependent on the location of the processing power of the server 52A-C. The processing power of the server 52A-C may be located anywhere within communication range and may be distributed in several locations. Cellular or landline telephones may be used to provide or to augment the radio signals 54 and/or 56.

FIG. 3A is a block diagram of the server 52A. The server 52A includes a radio transceiver 62, a virtual reference synthetic phase processor 63 and an anomaly detector 63A. The server 52A receives data for the reference positions (or it already has the reference positions) and the reference phases from the reference stations 51A-N. The processor 63 uses a virtual reference system (VRS) parametric model with the master synthetic vector 64 (in place of the virtual vector 27) together with the master and auxiliary reference network positions and phases for the reference network stations 51A-N and the three dimensional angles to the GPS satellites 16 for inferring the synthetic reference phases that would be measured if the measurements were made at the synthetic position 133 (instead of the virtual reference position 26).

The radio transceiver 62 transmits the radio signal 56 having synthesized reference data to the rover station 124A. The synthesized reference data still includes a correct geographical virtual reference position 26, as is conventional, but includes the synthetic reference phases in place of the actual or virtual reference phases that are used by a conventional rover station for positioning operation without the accuracy control of the present invention.

FIG. 3B is a block diagram of the server 52B. The server 52B includes the radio transceiver 62, a synthetic phase processor 65, a virtual reference processor 66 and an anomaly detector 66A. The server 52B receives data for the reference positions (or it already has the reference positions) and the reference phases from the reference stations 51A-N. The virtual reference processor 66 uses the virtual vector 27 and the reference positions and phases for the reference stations 51A-N with the three dimensional angles to the GPS satellites 16 for determining the virtual reference phases. The virtual reference processor 66 then passes the virtual reference phases and the virtual reference position 26 to the synthetic phase processor 65.

The synthetic phase processor 65 uses the synthetic offset vector 32 with the virtual reference phases and the three dimensional angles to the GPS satellites 16 for inferring the synthetic reference phases for the carrier phases that would be measured if the measurements were made at the synthetic position 133. The radio transceiver 62 transmits the radio signal 56 having synthesized reference data to the rover station 124B. The synthesized reference data still includes a correct geographical virtual reference position 26, as is conventional, but the phases for the GPS signals 14 are not the virtual reference phases that would be measured at the virtual reference position 26 but are instead synthetic reference phases that would be measured at the synthetic position 133.

FIG. 3C is a block diagram of the server 52C. The synthetic phase processor 65 is located separately from the virtual reference processor 66. The server 52C uses the public switch telephone network (PTSN) telephone system 68 for receiving reference data and uses the telephone system 68 for communication between the virtual reference processor 66 and the synthetic phase processor 65. The processor 65 receives data for the virtual reference position 26 and virtual reference phases from the virtual reference processor 66 and then infers the synthetic reference phases as described above. The processor 65 may be located adjacent to the rover station 124C having a local wired connection or a cellular telephone 69 may be used to pass the synthetic reference phases to the rover station 124C.

The virtual reference processors 63 and 66 determine double difference phase residuals between the master and auxiliary phases for current and previous phase measurements and pass the phase residuals to the respective anomaly detectors 63A and 66A. The anomaly detector 63A and 66A detects a phase residual anomaly when the phase residual is greater than a phase threshold corresponding to a selected distance or integrity limit 40 for RTK operation. The virtual reference processors 63 and the anomaly detectors 63A and 66A, respectively, may share hardware and software.

The rover GPS receivers 125A-C also determine double difference phase residuals. The phase residuals determined in the rover GPS receivers are the differences between the rover phases and the synthetic reference phases for current and previous phase measurements. The rover GPS receivers 125A-C pass the phase residuals to the respective anomaly detectors 126A-C. The anomaly detectors 126A-C also detect a phase residual anomaly when the phase residual is greater than a phase threshold corresponding to a selected distance or integrity limit 40 for RTK operation. The rover GPS receiver 125A-C and the anomaly detector 126A-C, respectively, may share hardware and software.

The integrity limit 40 corresponds to a zone about the rover position 38. When an anomaly is detected, the anomaly detector 63A, 66A or 126A-C inhibits the rover GPS receiver 125A-C from providing the rover position 38 to the user of the rover station 124A-C, or provides a notice to the user that the anomaly was detected and allows the user to decide whether or not to use the position 38. Alternatively, the rover station 124A-C provides a solution for the rover position 38 where the synthetic reference phase and measured rover phase for the particular GPS signal 14 associated with the anomaly are not used. The effect of the system 50 is that the rover position 38 has the controlled added position error 39 without degrading the integrity limit 40 of the RTK positioning solution for the rover position 38.

In the system 50, it may be beneficial to reduce the amount of data that is transmitted among various locations by sending differences between master and auxiliary reference network positions and/or reference network phases in place of the actual reference positions and phases. For example, the reference positions and phases for the auxiliary stations 51B-N may be transmitted as differences with respect to the reference position and phases the master reference station 51A.

FIG. 4 is a diagram showing a secure real time kinematic (RTK) GPS rover station 70 for receiving conventional reference data in a secure format from a GPS-based positioning system 71. The positioning system 71 includes at least one reference station 112 having a reference position that is established by a survey or some other means for receiving GPS signals 14 from GPS satellites 16. The reference station 112 measures the carrier phases of the GPS signals 14 and sends a radio signal 117 having secure reference data for the reference phases to the rover station 70. The security of the reference data may be maintained by the measures of the Digital Millennium Copyright Act of 1998 for preventing unauthorized access to a copyrighted work. Alternatively, the reference data may be encrypted.

The rover station 70 receives or generates or otherwise selects the synthetic offset vector 32. The synthetic offset vector 32 and the reference position of the reference station 112 define a synthetic position 33 as described above. The length and direction of the synthetic offset vector 32 is arbitrary but the length is normally a few meters or less. When it is desired for the rover station 70 to operate with existing RTK GPS-based reference systems, the reference station 112 may be a conventional reference station 12 described above with the addition of security measures for protecting the reference data from unauthorized access.

FIG. 4A is a block diagram of the rover station 70. The rover station 70 includes an RTK rover GPS receiver 74 including an anomaly detector 74A, a secure synthetic phase processor 75, and a radio transceiver 76. The rover GPS receiver 74 measures the carrier phases for the same GPS signals 14 that are measured by the reference station 112. The radio transceiver 76 may be replaced by a radio receiver without a transmitter if two-way communication is not required. A cellular telephone may be used for the radio transceiver 76. The radio transceiver 76 receives the reference position and the secure reference data for the reference phases in the radio signal 117.

The secure synthetic phase processor 75 selects the synthetic offset vector 32 and then uses the synthetic offset vector 32 with the reference position, the secure reference phases and the three dimensional angles to the GPS satellites for inferring the synthetic reference phases. The secure synthetic phase processor 75 performs processing on signals and data within physical boundaries of the processor 75 in a way that makes it difficult for an authorized user to alter the processing algorithms or view the signals or data. Further, the algorithms, signals, messages and data are protected by the access controls of the Digital Millennium Copyright Act of 1998.

The secure synthetic phase processor 75 passes the synthetic reference phases to the rover GPS receiver 74. The GPS receiver 74 uses the synthetic reference phases and the measured rover phases with the reference position and the three dimensional angles to the GPS satellites 16 to compute the rover position 38. The conventional rover station 18 would calculate the difference between the reference and rover phase measurements for the distance vectors represented by the vector d to the GPS satellite 16A. However, the rover station 70 arrives at estimated perpendicular distance vectors represented by d* instead of the vectors represented by d.

When the rover station 70 calculates its position with respect to the reference station 112, it arrives at a position 38 that has a vector position offset error 39 of equal length and in the opposite direction as the synthetic offset vector 32. The security measures in the secure processor 75 prevent the user from undoing the accuracy control of the present invention by using the measured reference phases instead of the synthetic reference phases. Measurements by the rover GPS receiver 74 from several GPS satellites 16 yield several perpendicular distance vectors d* and ultimately the position of the rover station 70 with respect to the reference station 112 with the added error 39. Using this technique the secure synthetic phase processor 75 is able to introduce the arbitrary added error 39 into the position 38 that is calculated by the rover station 70.

The rover GPS receiver 74 determines phase residuals from current and previous synthetic reference phases and measured rover phases and passes the phase residuals to the anomaly detector 74A. The anomaly detector 74A detects a phase residual anomaly when the phase residual is greater than a phase threshold corresponding to a selected integrity limit 40 for RTK operation. The integrity limit 40 corresponds to a zone about the rover position 38. When an anomaly is detected, the anomaly detector 74A inhibits the rover GPS receiver 74 from providing the rover position 38 to the user of the rover station 70, or provides a notification to the user that the anomaly was detected and allows the user to decide whether or not to use the position 38. Alternatively, the anomaly detector 74A provides a solution for the rover position 38 where the synthetic reference phase and measured rover phase for the particular GPS signal 14 associated with the anomaly are not used. The effect of the system 70 is that the rover position 38 has the controlled added position error 39 without degrading the integrity limit 40 of the RTK positioning solution for rover position 38.

FIG. 5 is a diagram showing a secure real time kinematic (RTK) GPS rover station 80A or 80B for receiving conventional reference data in a secure format from a network positioning system 81. The rover station 80A-B receives reference data having secure reference phases from the system 81 and selects the synthetic offset vector 32 for controlling the positioning accuracy that it provides. The security of the reference phases may be protected by the access control measures of the Digital Millennium Copyright Act of 1998 and/or by encryption. The positioning system 81 includes a server 123 and a network of reference network stations, referred to as 112A, 112B through 112N, having reference positions that are known from a survey or other means.

The reference network stations 112A-N include reference GPS receivers for receiving the GPS signals 14 from the GPS satellites 16 and measuring carrier phases. When it is desired for the rover station 80A-B to operate with existing RTK GPS-based reference systems, the reference stations 112A-N may be conventional reference stations 12A-N and the server 123 may be the conventional serve 23 described above with the addition of security measures for protecting the reference data. One of the reference stations, illustrated as 112A, may be designated as a master reference station and the other reference network stations 112B-N as auxiliaries.

The reference network stations 112A-N communicate with the server 123 in signals 122 and the server communicates with the rover station 80A-B with a radio signal 127 having a secure data format so that the reference phases cannot easily be used by an unauthorized user. The sum of the virtual vector 27 and the synthetic offset vector 32 is a master synthetic vector 64. The position of the master reference station 112A and the master synthetic vector 64 define a synthetic position 133.

FIG. 5A is a block diagram of the rover station 80A. The rover station 80A includes a radio transceiver 82, a secure virtual reference synthetic phase processor 83 and an RTK rover GPS receiver 84A including an anomaly detector 86A. The radio transceiver 82 receives data for the master and auxiliary reference positions and phases for the master and auxiliary reference network stations 112A-N in the radio signal 127. The radio transceiver 82 may be a radio receiver without a transmitter if two-way communication is not required. The radio transceiver 82 may be a cellular telephone. In order to reduce the amount of data that is transmitted, the reference positions and phases for the auxiliary reference stations 112B-N may be transmitted as differences from the reference position and phases of the master reference station 112A.

The synthetic phase processor 83 receives or generates or otherwise selects the synthetic offset vector 32 and then determines the virtual reference position 26, or negotiates with the server 123 to determine the virtual reference position 26. The virtual reference position 26 and the synthetic offset vector 32 define the synthetic position 133 where the synthetic position 133 is separated from the virtual reference position 26 by the synthetic offset vector 32. The processor 83 determines the master synthetic vector 64 from the vector sum of the virtual vector 27 and the synthetic offset vector 32 (or the virtual reference position 26 and the synthetic offset vector 32). The length and direction of the synthetic offset vector 32 are arbitrary but the length is normally a few meters or less.

The processor 83 then uses the master synthetic vector 64, in place of the virtual vector 27, together with the master and auxiliary reference network positions and phases for the reference network stations 112A-N and the three dimensional angles to the GPS satellites 16 for inferring the synthetic reference phases that would be measured if the measurements were made at the synthetic position 133. The processor 83 passes the synthetic reference data for the virtual reference position 26 and the synthetic reference phases to the rover GPS receiver 84A. The rover GPS receiver 84A measures the phases of the same GPS signals and uses the measured rover phases, the master and auxiliary reference positions and phases with the synthetic reference phases and the virtual reference position 26 for determining the rover position 38.

FIG. 5B is a block diagram of the rover station 80B. The rover station 80B is similar to the rover station 70 described above with the exception that the rover station 80B uses the virtual reference position 26 in place of the actual reference position of the reference station 112. The rover station 80B includes a radio transceiver 82, the rover GPS receiver 84B including an anomaly detector 86B, and a secure synthetic phase processor 85. The radio transceiver 82 receives data for the virtual reference position 26 and the virtual reference phases in the radio signal 127.

The processor 85 uses the synthetic offset vector 32 (or the difference between the virtual reference position 26 and the synthetic position 133) with the virtual reference position 26, virtual reference phases and the three dimensional angles to the GPS satellites for inferring the synthetic reference phases that would be measured at the synthetic position 133. The processor 85 passes the synthetic reference data for the virtual reference position 26 and the synthetic reference phases to the rover GPS receiver 84B. The rover GPS receiver 84B measures the phases of the same GPS signals and uses the measured rover phases with the synthetic reference phases and the virtual reference position 26 for determining the rover position 38.

The rover GPS receiver 84A-B determines phase residuals from current and previous synthetic reference phases and measured rover phases and passes the phase residuals to the anomaly detector 86A-B. The anomaly detector 86A-B detects a phase residual anomaly when the phase residual is greater than a phase threshold corresponding to a selected distance or integrity limit 40 for RTK operation. The integrity limit 40 corresponds to a zone about the rover position 38. When an anomaly is detected, the anomaly detector 86A-B inhibits the rover GPS receiver 84A-B from providing the rover position 38 to the user of the rover station 80A-B, or provides a notification to the user that the anomaly was detected and allows the user to decide whether or not to use the position 38. Alternatively, the anomaly detector 86A-B provides a solution for the rover position 38 where the synthetic reference phase and measured rover phase for the particular GPS signal 14 associated with the anomaly are not used.

The position 38 calculated by the rover station 80A-B relative to the virtual reference position 26 has the added positional offset error 39 that is equal in length and in the opposite direction as the synthetic offset vector 32. Using this technique the rover station 80A-B is able to arbitrarily introduce the added error 39 into the position 38 without degrading the integrity limit 40 of the RTK positioning solution for the rover position 38.

The secure synthetic phase processors 83 and 85 perform processing on signals and data that are embedded with the boundaries of the rover stations 80A and 80B in a way that makes it physically difficult for users of the rover stations 80A and 80B to alter the processing algorithms or view the signals or data. Further, the algorithms, signals, messages and data are protected by the access controls of the Digital Millennium Copyright Act of 1998.

FIGS. 6A and 6B are block diagrams of random reference generators referred to by reference numbers 90A and 90B, respectively. The random reference generator 90A is used in the synthetic phase processors 63 and 83 with a VRS parametric model for generating synthetic reference phases. The random reference generator 90B is used in the synthetic phase processors 35, 65, 75 and 85 with an actual or virtual reference position for generating synthetic reference phases. The synthetic reference phases are passed from the systems 30 and 50 to the RTK GPS receivers in the rover stations 118 or 124A-C for determining the rover position 38; or computed in secure processors within the rover station 70 or 80A-B for determining the rover position 38.

The random reference generators 90A and 90B include a random process vector generator 170. The random process vector generator 170 stores or receives values for a maximum rate of change and one or more maximum dimensions and uses the values as inputs to a random or pseudo-random process for continuously computing synthetic offset vectors 32. Importantly, because the synthetic offset vectors 32 are computed with a random or nearly random process, the added error 39 is not easily reversible by users or software programming in the rover stations.

The value or values for maximum dimensions may be a maximum radius for providing a spherical error zone, a maximum radius and a maximum length for providing a cylindrical error zone, three maximum lengths X, Y and Z for providing a box error zone, or the like. The error zones refer to a volume or a three dimensional range of the added error 39 for the rover position 38 about the position for the rover station that would be determined by an RTK rover station without accuracy control. For example, the added offset 39 for the box error zone has possible errors x, y and z in three dimensions of |x|≦X, |y|≦Y and |z|≦Z. The box error zone need not have equal or orthogonal dimensions. The values for the maximum dimensions z=0 or x and y=0 may be used to constrain the random process vector generator 170 so that the added error 39 is confined to horizontal or vertical directions, respectively.

The added error 39 may be of relatively large magnitude but low rate of change in any direction while the rover GPS receiver, constructed for fixed RTK operation, continues to use the resolved integer number of carrier phase cycles for its positioning. By continuing to use the integers, the rover position 38 has the integrity of the RTK GPS solution within the integrity limit 40 as small as a few centimeters even when the added error 39 is a few meters or more. The RTK rover position 38 has high integrity even when the accuracy is degraded with the present invention because the errors due to multipath are largely eliminated. It will be appreciated by those skilled in the art that merely dithering the reference carrier phase measurements directly and providing the dithered reference phases to the rover station may make it impossible for an RTK rover station to resolve the carrier phase integer, thereby losing the benefit of the high integrity of the RTK position solution.

The random reference generator 90A includes a vector summer 172 and a virtual reference phase synthesizer 174. The vector summer 172 sums the synthetic offset vector 32 with the virtual vector 27 for determining the master synthetic vector 64. The virtual reference phase synthesizer 174 uses the master synthetic vector 64 with the three dimensional angles to the GPS satellites 16, and the master and auxiliary reference positions and corresponding measured master and auxiliary reference carrier phases for the GPS signals 14 for computing the synthetic reference phases. The synthetic reference phases are then used as described above with the carrier phases measured by the RTK rover GPS receiver for computing the rover position 38.

The random reference generator 90B includes a phase synthesizer 175. The phase synthesizer 175 uses the synthetic offset vector 32 from the random process vector generator 170 with a reference carrier phase and the three dimensional angles to the GPS satellites 16 to compute the synthetic reference phases for the GPS signals 14. The reference carrier phase may be an actual reference phase measured at an actual reference position or a virtual reference phase computed for the virtual reference position 26. The synthetic reference phases are then used as described above with the carrier phases measured by the RTK rover GPS receiver for computing the rover position 38.

FIG. 7 is a system diagram showing a secure real time kinematic (RTK) GPS rover station 200A or 200B for operation with a reference station 212A or 212B, respectively, in a positioning system 201. The rover station 200A-B receives reference system data in a secure form from the system 201 and receives or generates a synthetic offset vector 232 with respect to the position of the reference station 212A-B. The secure reference data is used in the rover station 200A-B for computing a secure position 210. The rover station 200A-B then dithers the secure position 210 with the synthetic offset vector 232 to provide an unsecured rover position 238 having an added position error 239 to a user of the rover station 200A-B.

The vector for the added position error 239 is the same length and direction as the synthetic offset vector 232. Double difference phase residuals are monitored in the rover stations 200A-B for maintaining an integrity limit 240 for RTK operation about the secure position 210. The integrity limit 240 is an outer limit of a zone about the secure position 210. The added position error 239 offsets the secure position 210 to the unsecure position 238 without degrading the integrity limit 240 about the unsecure position 238. The integrity limit 240 may be twenty centimeters or less. The length and direction of the synthetic offset vector 232 and the added position error 239 are arbitrary but the length is normally a few meters or less.

The reference station 212A-B has a reference position that is established by a survey or some other means for receiving GPS signals 14 from GPS satellites 16 and measuring reference carrier phases. The synthetic offset vector 232 and the reference position define a synthetic position 233 of the reference station 212A-B. The reference station 212A-B sends a radio signal 217 having the secure reference data for the measured reference phases and the reference position to the rover station 200A-B. Information for the synthetic offset vector 232 may be included in the secure reference data transmitted to the rover station 200A or generated within the rover station 200B or received in the rover station 200A from some other secure source. When it is desired for the rover station 200B to operate with existing RTK GPS-based reference systems, the reference station 212B may be a conventional reference station 12 described above with the addition of security for the reference data that is transmitted to the rover station 200B.

FIG. 7A is a block diagram of an embodiment for the reference station 212A and the rover station 200A where the reference station 212A generates the synthetic offset vector 232. The reference station 212A includes a reference GPS receiver 252, a reference position memory 254, a synthetic vector generator 260, a secure reference data provider 262 and a radio transducer 264. The reference GPS receiver 252 receives and measures the carrier phases for the GPS signals 14. The reference position memory 254 stores the position of the reference station 212A. The synthetic vector generator 260 generates the synthetic offset vector 232. The secure data provider 262 processes the reference phases, reference position, and the synthetic offset vector 232 into a secure format for the reference data. The radio transceiver 264 issues the secure reference data in the radio signal 217 to the rover station 200A.

The rover station 200A includes a radio transceiver 272, an RTK rover GPS receiver 274 including an anomaly detector 275, and a position dither processor 277. The radio transceiver 272 may be replaced by a radio receiver without a transmitter if two-way communication is not required. A cellular telephone may be used for the radio transceiver 272.

The radio transceiver 272 receives the secure reference data for the synthetic offset vector 232 and a reference position and measured reference phases in the radio signal 217; and passes the reference position and phases to the GPS receiver 274 and the synthetic offset vector 232 to the position dither processor 277. The GPS receiver 274 measures the carrier phases for the GPS signals 14 for the same GPS satellites 16 and calculates the differences between the reference and rover phase measurements. The phase differences result in estimated perpendicular distance vectors represented by the vector d for the GPS signal 14A as described above for determining the secure rover position 210.

The position dither processor 277 dithers the secure rover position 210 with the synthetic offset vector 232 to provide the rover position 238 having the added position error 239. Preferably, the position dither processor 277 is a coded algorithm embedded in memory or signal processing hardware that is read or otherwise processed by the hardware and software in rover GPS receiver 274. Both the rover GPS receiver 274 and the position dither processor 277 must be secure from tampering by users of the rover station 200A in order to prevent users from undoing the accuracy control that is provided by the rover station 200A.

The rover GPS receiver 274 determines double difference phase residuals from current and previous reference phases and current and previous measured rover phases and passes the phase residuals to the anomaly detector 275. The anomaly detector 275 detects a phase residual anomaly when the phase residual is greater than a phase threshold corresponding to a selected distance for the integrity limit 240. The integrity zone about the secure position 210 is transferred by the position dither processor 277 to the integrity zone 240 about the dithered rover position 238. When an anomaly is detected, the anomaly detector 275 inhibits the position dither processor 277 from providing the rover position 238 to the user of the rover station 200, or provides a notification to the user that the anomaly was detected and allows the user to decide whether or not to use the position 238. Alternatively, the anomaly detector 275 provides a solution for the secure position 210 and the position dither processor 277 provides the rover position 238 where the reference phase and measured rover phase for the particular GPS signal 14 associated with the anomaly are not used.

FIG. 7B is a block diagram of an embodiment for the reference station 212B and the rover station 200B where the rover station 200B generates the synthetic offset vector 232. The reference station 212B and the rover station 200B operate as described above for the reference station 212A and rover station 200A with the exception that the synthetic offset vector 232 is generated by the synthetic vector generator 260 in the rover station 200B.

FIG. 8 is a system diagram showing a secure real time kinematic (RTK) GPS rover station 300A, 300B, 300C or 300D for operation with a server 323A, 323B, 323C or 323D, respectively, in a network positioning system shown generally with a reference identifier 301. The rover station 300A-D receives the secure reference system data from the system 301 and receives or generates the synthetic offset vector 232 for controlling the positioning accuracy that it provides. The synthetic offset vector 232 offsets the virtual reference position 26 for the system 301 to a synthetic position 333. The reference data is used by the rover station 300A-D for computing a secure position 310.

The rover station 300A-D then uses the synthetic offset vector 232 to dither the secure position 310 to provide an unsecure rover position 338 having the added position error 239 to a user of the rover station 300A-D. The integrity limit 240 for RTK operation represents the outer limit of a zone about the secure position 310. The added position error 239 offsets the secure position 310 to the unsecure position 338 without degrading the integrity limit 240 so that the integrity limit 240 becomes the outer limit of a zone about the position 338. The vector for the added position error 239 is the same length and direction as the synthetic offset vector 232. The length and direction of the synthetic offset vector 232 and the added position error 239 are arbitrary but the length is normally a few meters or less.

The positioning system 301 includes a network of reference network stations, referred to as 312A, 312B through 312N, having reference positions that are known from a survey or other means. The reference stations 312A-N measure the carrier phases of the GPS signals 14 from the GPS satellites 16 and send telephone or radio signals 322 having the reference system data for the measured phases to the server 323A-D. Information for the synthetic offset vector 232 may be included in the secure reference data transmitted to the rover station 300A,C or generated within the rover station 300B,D or received in the rover station 300A,C from some other secure source. The server 323A-D communicates with a radio signal 325 to send reference data in a secure format to the rover station 300A-D. One of the reference stations, illustrated as 312A, may be designated as a master reference station and the other reference network stations 312B-N as auxiliaries.

The system 301 determines the virtual reference position 26 and the virtual vector 27 from the master reference station 312A to the virtual reference position 26. When it is desired for the rover station 300B,D to operate with existing RTK GPS-based reference systems, the server 323B,D and the reference stations 312A-N may be the server 23 and the reference stations 12A-N of the prior art with the addition of the security for the reference data transmitted to the rover station 300B,D. It should be noted that the elements of the server 323A-D do not need to be in one physical location.

FIG. 8A is a block diagram of an embodiment for the server 323A and the rover station 300A where the server 323A generates the synthetic offset vector 232. The server 323A includes the synthetic vector generator 260, a VRS position phase processor 352 including an anomaly detector 353, a secure reference data provider 354, and a radio transceiver 356. The synthetic vector generator 260 generates and passes the synthetic offset vector 232 to the secure data provider 354.

The VRS reference position phase processor 352 receives the master and auxiliary reference phases from the reference stations 312A-N in the signal 322. The master and auxiliary reference positions are retained by the processor 352 or are received in the signal 322. The processor 352 uses the virtual reference position 26 with the master and auxiliary positions and phases from the reference stations 312A-N for determining virtual reference phases for the GPS signals 14 referred to the virtual reference position 26 and passes the reference data for the virtual reference position 26 and the virtual reference phases to the secure data provider 354. The anomaly detector 353 monitors double difference phase residuals between the current and previous measured phases between the reference stations 312A-N for preventing the reference data from being used by the rover station 300A when one of the current reference phases is outside a phase residual threshold corresponding to the integrity limit 240.

The secure data provider 354 processes the reference data into a secure format and passes the secure reference data to the radio transceiver 356. The radio transceiver 356 issues the secure reference data in the radio signal 325 to the rover station 300A. The rover station 300A includes the position dither processor 277, a radio transceiver 362, and a rover RTK GPS receiver 364 including an anomaly detector 365. The radio transceiver 362 receives secure reference data in the radio signal 325 and passes the reference system position and phase data to the rover GPS receiver 364 and the synthetic offset vector 232 to the position dither processor 277. The radio transceiver 362 may be a radio receiver without a transmitter if two-way communication is not required. The radio transceiver 362 may be a cellular telephone.

The rover GPS receiver 364 measures carrier phases for the GPS signals 14 for the same GPS satellites 16 as the reference GPS receiver 352 and then uses the reference system position and phase data for correcting the carrier phases that it measures and ultimately arrives at the secure position 310 with respect to the virtual reference position 26.

The rover GPS receiver 364 determines phase residuals from current and previous reference phases and measured rover phases and passes the phase residuals to the anomaly detector 365. The anomaly detector 365 detects a phase residual anomaly when the phase residual is greater than a phase threshold corresponding to a selected integrity limit 240. The integrity limit 240 corresponds to a zone about the rover position 310. The position dither processor 277 dithers the secure position 310 with the synthetic offset vector 232 for transferring the integrity limit 240 to the unsecure position 338 having the added error 239. When an anomaly is detected, the anomaly detector 365 inhibits the rover GPS receiver 364 from providing the rover secure position 310 to the position dither processor 277 and ultimately inhibits the rover station 300A from providing the rover position 338 to the user of the rover station 300A. Alternatively, the anomaly detector 365 provides the position 338 where the measured reference and rover phases for the particular GPS signal 14 associated with the anomaly are not used.

FIG. 8B is a block diagram of an embodiment for the server 323B and the rover station 300B where the rover station 300B generates the synthetic offset vector 232. The server 323B includes the VRS position phase processor 352 including the anomaly detector 353, the secure data provider 354, and the radio transceiver 356 as described above. The rover station 300B includes the synthetic vector generator 260, the position dither processor 277, the radio transceiver 362, and the rover RTK GPS receiver 364 including the anomaly detector 365 as described above.

The radio transceiver 356 transmits the reference system data for the virtual reference position 26 and the virtual reference phases in a secure format in the radio signal 325 to the rover station 300B. The synthetic vector generator 260 in the rover station 300B passes the synthetic offset vector 232 to the position dither processor 277. The position dither processor 277 dithers the secure position 310 with the synthetic offset vector 232 to provide the unsecure rover position 338 having the integrity limit 240 to the user of the rover station 300B.

FIG. 8C is a block diagram of an embodiment for the server 323C and the rover station 300C where the server 323C generates the synthetic offset vector 232. The server 323C includes the synthetic vector generator 260, a reference server processor 368, the secure data provider 354, and the radio transceiver 356. The synthetic vector generator 260 generates and passes the synthetic offset vector 232 to the secure data provider 354.

The reference server processor 368 receives the master and auxiliary reference phases from the reference stations 312A-N in the signal 322. The master and auxiliary reference positions are retained by the processor 368 or are received in the signal 322. The processor 368 passes the master and auxiliary reference positions and phases to the secure data provider 354.

The secure data provider 354 processes the reference data into a secure format and passes the secure reference data to the radio transceiver 356. The radio transceiver 356 issues the secure reference data in the radio signal 325 to the rover station 300C. The rover station 300C includes the position dither processor 277, the radio transceiver 362, and a rover RTK GPS receiver 374 including an anomaly detector 375. The radio transceiver 362 receives secure reference data in the radio signal 325 and passes the reference system position and phase data to the rover GPS receiver 374 and the synthetic offset vector 232 to the position dither processor 277. The radio transceiver 362 may be a radio receiver without a transmitter if two-way communication is not required. The radio transceiver 362 may be a cellular telephone.

The rover GPS receiver 374 measures carrier phases for the GPS signals 14 for the same GPS satellites 16 as the reference stations 312A-N and then uses the reference system position and phase data for correcting the carrier phases that it measures and ultimate arrives at the secure position 310 with respect to the virtual reference position 26. The secure position 310 is passed to the position dither processor 277.

The rover GPS receiver 374 determines phase residuals from current and previous reference phases and measured rover phases and passes the phase residuals to the anomaly detector 375. The anomaly detector 375 detects a phase residual anomaly when the phase residual is greater than a phase threshold corresponding to a selected integrity limit 240. The integrity limit 240 corresponds to the outer limit of a zone about the rover position 310. The position dither processor 277 dithers the secure position 310 with the synthetic offset vector 232 for transferring the integrity limit 240 to the unsecure rover position 338 having the added position error 239 with respect to the virtual reference position 26. When an anomaly is detected, the anomaly detector 375 inhibits the rover GPS receiver 374 from providing the rover secure position 310 to the position dither processor 277 and ultimately inhibits the rover station 300C from providing the rover position 338 to the user of the rover station 300C. Alternatively, the anomaly detector 375 provides a solution for the rover position 338 where the reference phase and measured rover phase for the particular GPS signal 14 associated with the anomaly are not used.

FIG. 8D is a block diagram of an embodiment for the server 323D and the rover station 300D where the rover station 300D generates the synthetic offset vector 232. The server 323D includes the reference server processor 368, the secure data provider 354, and the radio transceiver 356 as described above. The rover station 300D includes the synthetic vector generator 260, the position dither processor 277, the radio transceiver 362, and the rover RTK GPS receiver 374 including the anomaly detector 375, as described above.

The radio transceiver 356 transmits the reference system data for the master and auxiliary reference positions and phases (or difference between the master and auxiliary positions and phases) in a secure format in the radio signal 325 to the rover station 300D. The synthetic vector generator 260 in the rover station 300D generates and passes the synthetic offset vector 232 to the position dither processor 277. The position dither processor 277 dithers the secure position 310 with the synthetic offset vector 232 to provide the unsecure rover position 338 to the user of the rover station 300D.

Using these techniques the secure rover station 200A-B,300A-D is able to introduce the arbitrary added error 239 into the position 238,338 that is provided by the rover station 200A-B,300A-D without degrading the integrity limit 240 that is computed for the secure position 210,310. The reference data must be secure and both the rover GPS receiver 274,364,374 and the position dither processor 277 must be secure from tampering by users in order to prevent users from undoing the accuracy control that is provided by the rover station 200A-B,300A-D.

Anomaly detectors, such as 63A, 66A, 74A, 86A-B, 120, 126A-C, 275, 353, 365 and 375, are described above for detecting anomalies (also known as is outliers) for phase residuals on a satellite-by-satellite basis for real time kinematic (RTK) position determinations. The anomaly is detected when the phase residual exceeds a selected phase residual limit. The phase residual limit is selected so that when the phase residuals are within the phase residual limit, the position determination has the designated integrity limit 40,240.

The algorithms, signals, messages and data in the rover GPS receiver 274,364,374, the position dither processor 277, and the reference data in the signals 127, 217 and 325 are provided with the access control measures of the Digital Millennium Copyright Act of 1998. The reference data in the signals 127, 217 and 325 may also be protected by encryption. The rovers 200A-B and 300A-D perform processing on signals and data that are embedded with the boundaries of the rovers 200A-B and 300A-D in a way that makes it mechanically or electrically difficult for an unauthorized user of the rovers 200A-B and 300A-D to alter the algorithms or view the signals or data. All users are unauthorized users unless they are designated by the provider of the algorithms, signals or data as authorized users.

It may be noted that the positioning system 30, 50, 71, 81, 201 or 301 could be used as the basis for a fee-based RTK GPS service wherein the price for the service is based upon the accuracy of the positioning.

FIG. 9 is a block diagram illustrating the synthetic vector generator 260 and the position dither processor 277. The synthetic vector generator 260 includes a random process vector generator 380. The random process vector generator 380 stores or receives values for a maximum rate of change and one or more maximum dimensions and uses the values as inputs to a random or pseudo-random process for continuously computing the synthetic offset vector 232. The position dither processor 277 includes a summer 382 for summing the synthetic offset vector 232 with the secure position 210, 310. Importantly, because the synthetic offset vector 232 is the same as the added positional error 239 and the synthetic offset vector 232 is computed with a random or nearly random process, the added positional error 239 is not easily reversible by unauthorized users.

The value or values for maximum dimensions may be a maximum radius value for providing a spherical error zone, a maximum radius and a maximum length for providing a cylindrical error zone, three maximum lengths X, Y and Z for providing a box error zone, or the like. The error zones refer to a volume or a three dimensional range of the added position error 239 for the dithered (unsecure) rover position 238,338 about the secure rover position 210,310. For example, the added error positional error 239 for the box error zone has possible errors x, y and z in three dimensions of |x|≦X, |y|≦Y and |z|≦Z. The box error zone need not have equal or orthogonal dimensions. The values for the maximum dimensions z=0 or x and y=0 may be used to constrain the random process vector generator 380 so that the added position error 239 is confined to horizontal or vertical directions, respectively.

The added error 239 may be of relatively large magnitude in any direction while the rover GPS receiver, constructed for fixed RTK operation, continues to use the resolved integer number of carrier phase cycles for its positioning. By continuing to resolve the integers, the rover position 238,338 has the integrity of the RTK GPS solution within the integrity limit 240 as small as a few centimeters even when the added error 239 is a few meters or more. The RTK rover position 238,338 has high integrity when the accuracy is degraded because the errors due to multipath are largely eliminated. It will be appreciated by those skilled in the art that merely dithering the reference position directly and providing the dithered reference position to the rover could make it impossible for an RTK rover station to resolve the carrier phase integer, thereby losing the benefit of the high integrity of the RTK position solution.

FIG. 10 is a flow chart of steps of a method for providing the synthetic reference phases from the reference system 30 having one reference station to one or more rover stations. The steps may be embodied in a tangible medium 600 containing instructions that may be read by a processor or processors for causing the system to carry out the steps. The medium 600 may be constructed with one or more memory devices such as compact disks, electronic memory chips, hard disks, digital video devices, or the like. The processor may be a device commonly known as a computer or a microprocessor.

A synthetic offset vector is received or generated or otherwise selected in a step 602. In a step 604, GPS signals are received at the reference station by a real time kinematic (RTK) GPS receiver. In a step 606 the reference GPS receiver measures carrier phases of the GPS signals at a reference position.

A synthetic position is defined by the reference position and the synthetic offset vector. In a step 608 the reference system uses the synthetic offset vector and the measured reference phases for determining synthetic reference phases for the GPS signals that would be received at the synthetic position. In a step 612 the reference system transmits synthesized reference data that includes the synthetic reference phases to the rover station.

A GPS rover station having an RTK GPS receiver receives the synthesized reference data in a step 614. In a step 616 the rover GPS receiver receives the GPS signals from the same GPS satellites. In a step 618 the rover GPS receiver measures the carrier phases of the same GPS signals. In a step 622, the synthetic reference phases and the measured rover phases are used for testing the integrity of the phase measurements. In a step 624, when integrity has been verified, the rover station uses the reference position, the synthetic reference phases and the rover phases for determining its position. The position that is determined by the rover station has the same RTK integrity as if it were determined with the phases for the reference position but with an added offset error, unknown to the rover station, equal in length to the synthetic offset vector.

FIG. 11 is a flow chart of steps of a method for providing the synthetic reference phases from a reference network system, such as the reference network system 50, to one or more rover stations. The steps may be embodied in a tangible medium 650 containing instructions that may be read by a processor or processors for causing the system to carry out the steps. The medium 650 may be constructed with one or more memory devices such as compact disks, electronic memory chips, hard disks, digital video devices, or the like. The processor may be a device commonly known as a computer or a microprocessor.

A synthetic offset vector is received or generated or otherwise selected in the step 602. In a step 652 one of the network of reference stations is designated as the master reference station. In a step 654 GPS signals are received at the reference network stations. In a step 656 the real time kinematic (RTK) GPS receivers at the reference stations measure reference network phases for the carriers of the GPS signals at reference network positions.

A virtual reference position is selected by the system or negotiated between the system and the rover station in a step 662. In a step 664 a virtual vector is calculated between the position of the master reference station and the virtual reference position. In a step 666 a master synthetic vector is calculated by adding the virtual vector and the synthetic offset vector. A synthetic position is defined by the virtual reference position and the synthetic offset vector or equivalently by the position of the master reference station and the master synthetic vector. In a step 674 the system uses the master synthetic vector, the reference network positions and the measured reference network phases for mathematically determining the synthetic reference phases that would be measured for GPS signals received at the synthetic position. In a step 675 the reference system uses double difference phase residuals of the master and auxiliary reference phases for testing integrity. In a step 676 when the integrity of the reference phases has been determined, the system transmits synthesized reference data including the synthetic reference phases. In a step 678 the rover station receives the synthetic reference data.

A GPS rover station having a real time kinematic (RTK) GPS receiver receives the GPS signals from the same GPS satellites in a step 682. In a step 684 the rover GPS receiver measures the carrier phases of the GPS signals. In a step 685, the synthetic reference phases and the measured rover phases are used for testing the integrity of the phase measurements. In a step 686, when integrity has been verified, the rover station uses the synthetic reference phases for determining its position. The accuracy of the position that is determined by the rover station has the same integrity as if it were determined in an RTK solution with the phases for virtual reference position but with an added offset error, unknown to the rover station, equal in length to the synthetic offset vector.

FIG. 12 is a flow chart of steps of a method for computing synthetic reference phases in a rover station, such as the rover station 70, and then using the synthetic reference phases for computing a rover position having an added error. The steps may be embodied in a tangible medium 700 containing instructions that may be read by a processor or processors for causing the rover station to carry out the steps. The medium 700 may be constructed with one or more memory devices such as compact disks, electronic memory chips, hard disks, digital video devices, or the like. The processor may be a device commonly known as a computer or a microprocessor.

A synthetic offset vector is received or generated or otherwise selected in a secure synthetic phase processor in the rover station in a step 702. In a step 704, GPS signals are received at a reference position by a real time kinematic (RTK) GPS receiver in a GPS reference station. A synthetic position is defined by the reference position and the synthetic offset vector. In a step 706 the reference GPS receiver measures carrier phases of the GPS signals. In a step 708 the reference system transmits a signal having secure reference data that includes the reference position and measured reference phases.

The rover GPS receiver receives the secure reference data from the system in a step 714. In a step 716 the rover GPS receiver receives the GPS signals from the same GPS satellites. In a step 718 the rover GPS receiver measures the carrier phases of the GPS signals. In a step 722 the rover station uses the reference position and phases and the synthetic offset vector for inferring the synthetic reference phases for the GPS signals that would be received at the synthetic position. In a step 723, the synthetic reference phases and the measured rover phases are used for testing the integrity of the phase measurements. In a step 724, when integrity has been verified, the rover GPS receiver uses the reference position, the synthetic reference phases and the measured rover phases for determining its position. The accuracy of the position that is determined by the rover station has the same RTK integrity as if it were determined with the true reference phases but with an added offset error, unknown to the rover station, equal in length to the synthetic offset vector.

FIG. 13 is a flow chart of steps of a method for computing synthetic reference phases in a rover station such as the rover station 80A or 80B, and then using the synthetic reference phases for computing a rover position having an added error. The steps may be embodied in a tangible medium 750 containing instructions that may be read by a processor or processors for causing the rover station to carry out the steps. The medium 750 may be constructed with one or more memory devices such as compact disks, electronic memory chips, hard disks, such as digital video devices, or the like. The processor may be a device commonly known as a computer or a microprocessor.

A synthetic offset vector is received or generated or otherwise selected in a secure phase processor in the rover station in the step 702. In a step 752 one of a network of reference stations is selected or designated as the master reference station. In a step 754 GPS signals are received by real time kinematic (RTK) GPS receivers at GPS reference network stations. In a step 756 the reference GPS receivers measure reference network phases for the carriers of the GPS signals.

A virtual reference position is selected by the system or negotiated between the system and the rover station in a step 762. In a step 764 a virtual vector is calculated from the position of the master reference station to the virtual reference position. A synthetic position is defined by the virtual reference position and the synthetic offset vector. In a step 774 the system or the rover station computes virtual reference phases from the virtual reference vector and the reference network positions and phases. In a step 776 the system transmits secure reference data having the measured master and auxiliary reference positions and phases or the virtual reference position and phases to the rover station.

The rover GPS receiver receives the GPS signals from the same GPS satellites in a step 778. In a step 782 the rover GPS receiver measures the carrier phases of the GPS signals. In a step 784 the rover station uses the synthetic offset vector directly or indirectly through the master synthetic vector (calculated by adding the virtual vector and the synthetic offset vector), and the virtual reference position and phases or the master and auxiliary positions and phases for computing the synthetic reference phases for the synthetic position. In a step 785, the synthetic reference phases and the measured rover phases are used for testing the integrity of the phase measurements. In a step 786, when integrity has been verified, the rover GPS receiver uses the virtual reference position and phases and the synthetic reference phases for determining a position. The accuracy of the position that is determined by the rover station has the same RTK integrity as if it were determined for the virtual reference phases but with an added offset error, unknown to the rover station, equal in length to the synthetic offset vector.

FIG. 14 is a flow chart of steps of a method for adding an error, shown and described above as the added error 239, to a rover position. The steps may be embodied in a tangible medium 800 containing instructions that may be read by a processor or processors for causing the rover station to carry out the steps. The medium 800 may be constructed with one or more memory devices such as compact disks, electronic memory chips, hard disks, digital video devices, or the like. The processor may be a device commonly known as a computer or a microprocessor.

A synthetic offset vector is received or generated or otherwise selected in the rover station in a step 802. In a step 804, GPS signals are received at a reference position by a real time kinematic (RTK) GPS receiver in a GPS reference station. In a step 806 the reference GPS receiver measures carrier phases of the GPS signals. In a step 808 the reference system transmits secure reference data that includes the reference position and measured reference phases in a secure format to the rover station.

The secure rover station receives the reference data in a step 814. The rover GPS receiver receives the GPS signals from the same GPS satellites in a step 816. In a step 818 the rover GPS receiver measures the carrier phases of the same GPS signals. In a step 820, the reference and rover phases are tested for the integrity of the phase measurements. In a step 822 the rover GPS receiver uses the reference data and the rover GPS phase measurements for determining a secure rover position. The secure position is not made available to normal users of the rover station. In a step 824 the rover station dithers the secure rover position with the synthetic offset vector from the step 802 to provide an unsecure rover position to users of the secure rover station. The unsecure rover position has the added error equivalent to the synthetic position offset.

FIG. 15 is a flow chart of steps of a method for adding an error, shown and described above as the added error 239, to a rover position. The steps may be embodied in a tangible medium 850 containing instructions that may be read by a processor or processors for causing the rover station to carry out the steps. The medium 850 may be constructed with one or more memory devices such as compact disks, electronic memory chips, hard disks, digital video devices, or the like. The processor may be a device commonly known as a computer or a microprocessor.

A synthetic offset vector is received or generated or otherwise selected in a secure rover station in the step 802. In a step 852 one of a network of reference stations is selected or designated as the master reference station. In a step 854 GPS signals are received by real time kinematic (RTK) GPS receivers at GPS reference network stations. In a step 856 the reference GPS receivers measure reference network phases for the carriers of the GPS signals.

A virtual reference position is selected by the system or negotiated between the system and the rover station in a step 862. In a step 864 a virtual vector is calculated from the position of the master reference station to the virtual reference position. In a step 874 the system or rover computes virtual reference phases from the virtual reference vector and the reference network positions and phases. In a step 876 the system transmits reference data including the reference phases in a secure format to the rover station. In a step 877 the secure rover station receives the reference data.

A GPS receiver in the rover station receives the GPS signals from the same GPS satellites in a step 878. In a step 882 the rover GPS receiver measures the carrier phases of the same GPS signals. In a step 883, the reference network and rover phases are tested for the integrity of the phase measurements. In a step 884 the rover GPS receiver uses the reference network and the rover GPS phase measurements for determining a secure rover position. The secure position is not made available to normal users of the rover station. In a step 886 the rover station dithers the secure rover position with the synthetic offset vector to provide an unsecure rover position to users of the rover station. The unsecure rover position has the added error that is equivalent to the synthetic offset vector.

The reference data, rover stations, and structural parts of the rover stations are described in several embodiments of the invention as secure. In the context of the invention, the term “secure” means that security provisions have been made to inhibit or prevent unauthorized users from viewing, accessing or changing the signals, the data or the algorithms in the secure elements. The security provisions may include encryption, the privacy measures of the Digital Millennium Copyright Act of 1998 for preventing unauthorized access to a copyrighted work, and physical constraints such as sealed packaging and using embedded codes, signals and data in ways so that it is physically or mechanically difficult view, access or change the codes, signals and data. A secure position is only available to an authorized user. An unsecure position is a position that is available to unauthorized (normal) users.

The provider of the reference system, the provider of the reference data, and/or the provider of the algorithms for controlling the positional accuracy of the rover station designate the parties that are authorized users. All other users of the reference system, reference data and/or rover stations are unauthorized users. Anyone not designated by the provider for controlling the position accuracy is an “unauthorized user”. The unauthorized user is constrained from using the rover station for obtaining a positional accuracy that does not have the added positional error that is controlled by the provider. In general the provider is or represents the seller of the reference system, the reference data, or the rover station and the unauthorized user is the one who is using the rover station for field work as the normal user or end user for the rover positions.

The FIGS. 2 and 10 illustrate embodiments where a single reference station system generates synthetic reference phases and provides the synthetic reference phases in an unsecure (public) signal to an existing rover station. The FIGS. 3 and 11 illustrate embodiments where a reference network system generates synthetic reference phases and provides the synthetic reference phases in a public signal to an existing rover station. The FIGS. 4 and 12 illustrate embodiments where a secure rover station synthesizes synthetic reference phases from the true reference phases received in a secure (private) signal from a single reference station system. The FIGS. 5 and 13 illustrate embodiments where a secure rover station synthesizes synthetic reference phases from the true reference phases received in a secure (private) signal from a reference network system.

The FIGS. 7 and 14 illustrate embodiments where a secure rover station receives reference data in a secure (private) signal from a single reference station system, computes a secure private true position and dithers the secure position with a synthetic offset vector for providing an unsecure position to a user. The FIGS. 8 and 15 illustrate embodiments where a rover station receives reference data in a secure (private) signal from a reference network system, computes a true secure private position and dithers the secure position with a synthetic offset vector for providing an unsecure position to a user.

FIG. 16 is a block diagram of a global navigation satellite system (GNSS) system 950 having real time kinematic (RTK) reference system 960, a rover receiver 970, and a computer apparatus 1000. The computer apparatus 1000 receives the GNSS RTK reference data including GNSS satellite reference carrier phases from the reference system 960 through a communication path 965; and receives the GNSS RTK rover data including GNSS satellite rover carrier phases from the rover receiver 970 through a communication path 975.

The path 965 may be the internet, a radio signal, or some other signal medium or a tangible medium. The path 975 may be a tangible medium such as a flash card, a docking of the rover receiver 970 to the computer apparatus 1000, or some other tangible medium or a signal medium. The computer apparatus 1000 has a post-processing capability for determining a position 1010 (FIG. 17) of the rover receiver 970 that is kept secure within the apparatus 1000 and a position 1038 (FIG. 17) that is made available to a user of the apparatus 1000 where the user-available position 1038 has a selectable degraded accuracy or precision 1041 (FIG. 17).

In general the computer apparatus 1000 is located in a different place than the rover receiver 970 and post processes the reference and rover data at a later time than the GNSS signals are received for the determinations of the reference and rover carrier phases. The geographical positions of the rover receiver 970 and the computer apparatus 1000 are entirely unrelated. The later time for processing the data is typically hours or days but possibly years. In a typical operation the rover receiver 970 has been moved to another site before the apparatus 1000 determines the secure rover position 1010. In an embodiment the computer apparatus 1000 is greater than 150 kilometers from the rover receiver 970 and the post processing is more than 12 hours after the determinations of the reference and rover carrier phases.

FIG. 17 is a two dimension geographical position diagram representing a three dimension positions of the rover receiver 970 with respect to the reference system 960. A secure rover position 1010 is determined by a computer apparatus 1000 having a centimeter level kinematic post-processing (KPP) capability. The secure position 1010 has an intrinsic uncertainty referred to as an intrinsic accuracy or an intrinsic precision 1040 due to the GNSS satellites that are used to compute the secure position 1010 and several other factors such as signal to noise ratios and/or estimated multipath errors of the satellite signals. The secure position 1010 is analogous to the secure positions 210 (FIGS. 7, 7A-B and 9) and 310 (FIGS. 8, 8A-D and 9) and the intrinsic precision 1040 is analogous to the integrity limits 40 (FIGS. 2-5) and 240 (FIGS. 7-8) described above.

The user-available rover position 1038 is provided to a user with an added error 1039 that degrades the intrinsic precision 1040. The stream of added errors 1039 have a standard deviation to provide the user-available rover position 1038 at the selected precision 1041. The added errors 1039 are in general 3 dimensional vectors having varying directions and lengths. The user-available rover position 1038 is analogous to the rover positions 38 (FIGS. 2, 3, 3A-C, 4, 4A, 5 and 5A-B), 238 (FIGS. 7, 7A-B and 9), and 338 (FIGS. 8, 8A-D and 9) and the added error 1039 is analogous to added errors 39 (FIGS. 2-5) and 239 (FIGS. 7-8) described above.

The intrinsic precision 1040 and the selected precision 1041 are shown as circles having diameters representative of standard deviations. In the illustrated case, an instantaneous user-available rover position 1038 is within the standard deviation of the selected precision 1041. The added error 1039 may sometimes place the user rover position 1038 outside the standard deviation of the selected precision 1041.

FIG. 18 is a block diagram of the computer apparatus 1000 having a specialized kinematic post-processing (KPP) capability. The computer apparatus 1000 may be a standard personal computer with the addition of software and/or hardware for a kinematic post-processor (KPP) 1102, a vector offset generator 1104 and a position dither processor 1106. The post-processor 1102 receives RTK reference data including measurements or determinations of GNSS reference carrier phases from the GNSS RTK reference system 960 and RTK rover data including measurements or determinations of GNSS rover carrier phases from the RTK rover receiver 970.

The post-processor 1102 uses the RTK reference and rover data to compute the secure rover position 1010. The vector offset generator 1104 generates offset vectors for providing the added errors 1039. The position dither processor 1106 adds the offset vectors to the secure position 1010 to provide the user-available rover position 1038 having the selected precision 1041.

The post-processor 1102 includes a dilution of precision (DOP) calculator 1108 and an intrinsic precision estimator 1109. The DOP calculator 1108 calculates the horizontal DOP (HDOP) and the vertical DOP (VDOP) for the geometry of the satellites that are used in the calculation of the secure position 1010. The intrinsic precision estimator 1109 estimates the intrinsic precision 1040 of the secure position 1010. In an embodiment the intrinsic precision 1040 is computed with a covariance matrix of the position fix for the secure position 1010.

The vector offset generator 1104 includes an accuracy leveler 1112 and a random process generator 1114. The accuracy leveler 1112 receives information for the selected precision 1041 and uses the intrinsic precision 1040 to compensate the selected precision 1041 to compute a dither level. The compensation reduces or eliminates variation of the selected precision 1041 due to variation of the intrinsic precision 1040. The effect of the compensation is to hold the selected precision 1041 approximately constant and independent of the intrinsic precision 1040 as long as the intrinsic precision 1040 is better than the selected precision 1041.

The random process generator 1114 generates a sequence for a stream of offset vectors having a standard deviation of the dither level that was computed by the accuracy leveler 1112. In an embodiment the random process generator 1114 has east, north and up random processors 1114E, 1114N and 1114U, respectively, for providing the east, north and up (vertical) components of the offset vectors. In an embodiment the random process generator 1114 uses a selected settling time constant, a selected bias, and a seed for generating the sequence of offset vectors. The position dither processor 1106 adds the offset vectors to the secure rover position 1010 to provide the user-available rover position 1038 with the selected precision 1041 to the user.

FIG. 19 is a block diagram of the random process generator 1114 having a seed generator 1120, a bias generator 1121, a sequence generator 1122 and an elongator 1124. The seed generator 1120 generates a seed value to be associated with a particular set of GNSS RTK reference and rover data that is received by the computer apparatus 1000 for a particular KPP position fix. The bias generator 1121 generates a position error bias to be associated with the particular data set.

The sequence generator 1122 uses the seed value to generate or select a first pseudorandom sequence w(i) having a near Gaussian distribution. In an embodiment the white noise w(i) sequence can be generated by first taking a uniformly distributed random number (C function-rand ( )), converting it to the range 0-1, then normalizing it via an approximation to the Error Function (see Abramovitz and Stegan, 1970, Handbook of Mathematical Functions, 9th edition, Dover, N.Y., for approximations to the Error Function). This process is repeatable by seeding the rand( ) function as this function is pseudorandom and generates the variable using a known algorithm. Using the same seed each time results in the same sequence of uniform value which in turn results in the same sequence of w(i). Finally, using the same w(i) then results in the same g(i). In an embodiment the seed value is effectively the first 32 bits as stored in the header of the file for GNSS reference and rover data set.

By always using the same seed value, the sequence w(i) is repeatable so that post-processing the same data set multiple times gives the same user-available rover position 1038. This is beneficial because otherwise a user would have possibly disconcerting results of seeing different position for the position 1038 each time the post processing was performed.

The elongator 1124 uses a selected time constant T with the sequence w(i) and the dither level σ from the accuracy leveler 1112 and the selected bias error k₀ to generate a second (elongated) pseudorandom sequence g(i) for the offset vectors. The elongated sequence g(i) has the selected time constant T to provide a longer time to average to the bias error k₀, or to average to zero if the bias is zero, than the first sequence w(i). The elongator 1124 passes the offset vectors to the position dither processor 1106. The use of the mean bias error is beneficial so that the user-available position 1038 is not averaged to the secure position 1010 by running the post processing for a longer time.

Elongation Process

The process for elongating a pseudorandom sequence for degrading the secure position 1010, which may be a KPP fixed solution, is via position coordinate dithering, i.e. applying pseudo random errors to the latitude, longitude and height components of the rover position 1010. A first-order Gauss-Markov (GM) process may be used because the GM process produces a distribution of errors that follow a normal distribution when averaged over time. The GM process can be controlled so that the errors change slowly over time and therefore the dithering process would be difficult to remove by the user as opposed to pure Gaussian errors of the first sequence w(i) which could be quickly averaged out.

An algorithm 1 shows a GM process for generating the elongated sequence for the stream of offset vectors.

$\begin{matrix} {{g(i)} = {{^{(\frac{{- \Delta}\; t}{T})}{g\left( {i - 1} \right)}} + {{w(i)}\sigma \sqrt{\left( {1 - ^{(\frac{{- 2}\; \Delta \; t}{T})}} \right)}}}} & (1) \end{matrix}$

Where:

g(i) is the value of the GM process at epoch i. Δt time difference between samples. T correlation time of the GM process. w(i) white noise at epoch i which has a Normal distribution. e Euler's number. σ Dither level standard deviation.

The selected correlation time constant T of the GM process defines how quickly the dithering errors can be averaged out, while σ defines the magnitude of the GM dithering errors. The GM dithering errors would take roughly 3 times the time constant T to average out.

The white noise w(i) input can be generated by first taking a uniformly distributed random number (C function-rand( )), converting it to the range 0-1, is then normalizing it via an approximation. The dither process g(i) mean of the algorithm 1 approaches zero average for long time periods (more than about 3T).

A post processing user could conceivably negate this dither process with static data by processing his data set but would need a sufficiently long time more than about 3T to do so.

An algorithm 1a shows a position error bias k₀ for generating a random process dither sequence b₀(i) that may be used for providing the sequence of offset vectors.

b ₀(i)=k ₀ +w(i)σ  (1a)

Where:

b₀(i) is the value of the output process at epoch i. w(i) white noise at epoch i which has a Normal distribution. k₀ is a bias of the dither sequence g(i). σ Dither level standard deviation.

The dither process b₀(i) mean never averages to zero so it cannot be negated for any length of time. The bias k₀ may be a variable determined randomly for each data set. To avoid confusing a user the k₀ might be limited to small changes for small user edits of the data set.

An algorithm 1b shows a combination of the ideas of the GM process for generating an elongated sequence g(i) and the position bias error k₀ for generating the stream of offset vectors.

$\begin{matrix} {{g(i)} = {{^{(\frac{{- \Delta}\; t}{T})}{g\left( {i - 1} \right)}} + {k_{0}\left\lbrack {1 - ^{({1 - \frac{\Delta \; t}{T}})}} \right\rbrack} + {{w(i)}\sigma \sqrt{\left( {1 - ^{(\frac{{- 2}\Delta \; t}{T})}} \right)}}}} & \left( {1b} \right) \end{matrix}$

Where:

g(i) is the value of the GM process at epoch i. Δt time difference between samples. T correlation time of the GM process. k₀ is a bias of the dither sequence g₀(i). w(i) white noise at epoch i which has a Normal distribution. e Euler's number. σ Dither level standard deviation.

The algorithm 1b shows a pseudorandom sequence g(i) computed with the bias error k₀ to provide a non zero average with a standard deviation σ for the computed dither level. The dither process g(i) mean for the algorithm 1b approaches the bias k₀ for long time periods more than about 3T.

Standard Deviation Calculation for Dither

FIG. 20 is a block diagram of an embodiment of the position dither processor 1106 and the vector offset generator 1104 where the vector offset generator 1104 includes the accuracy leveler 1112 and the elongator 1124. In this embodiment, covariance information for the secure position 1010 is used to compensate the selected precision 1041 with the intrinsic precision 1040 in order to decrease the lengths of the offset vectors as the intrinsic precision 1040 gets worse and increase the lengths of the offset vectors as the intrinsic precision 1040 gets better so that the selected precision 1041 that is provided to the post processing user remains about the same even when the intrinsic precision 1040 varies.

The accuracy leveler 1112 includes a horizontal precision compensator 1132, a vertical precision compensator 1134, an east horizontal scaler 1136, a north horizontal scaler 1138, and east, north and up square root functions 1139, 1140 and 1141, respectively. The selected precision 1041 may be composed of a selected horizontal precision (hzPrecInput) and a selected vertical precision (vtPrecInput). The horizontal and vertical precisions may be selected independently. The accuracy leveler 1112 receives information for the selected horizontal precision and the selected vertical precision respectively, and issues dither levels for east, north and up.

The horizontal precision compensator 1132 compensates the level of the hzPrecInput with east and north covariance diagonal components Qee and Qnn, respectively, to provide a compensated horizontal dither level variance (hzVarDither). Similarly, the vertical precision compensator 1134 compensates the vtPrecInput with the up covariance diagonal component Quu to provide a compensated vertical dither level variance (vtVarDither).

The east horizontal scaler 1136 uses the covariance components Qee and Qnn to scale the hzVarDither to provide an east dither level variance (EastVarDither). Similarly, the north horizontal scaler 1138 uses the covariance components Qee and Qnn to scale the hzVarDither to provide a north dither level variance (NorthVarDither). The east, north and up square root functions 1139, 1140 and 1141 compute the square roots of EastVarDither, NorthVarDither and vtVarDither, respectively, to provide dither level standard deviations for east σ_(e), north σ_(n), and up σ_(u), respectively.

The elongator 1124 includes an east sequence elongator 1144E, a north sequence elongator 1144N, and a vertical sequence elongator 1144U. In an embodiment the epochs i of the sequence w(i) have epoch offsets to provide an east sequence w(i_(e)), a north sequence w(i_(n)) and an up sequence w(i_(u)). Alternatively, three separate sequences w_(e)(i), w_(n)(i), and w_(u)(i) may be used for east, north and up, respectively. The east sequence elongator 1144E processes the east dither level standard deviation σ_(e) with the sequence w(i_(e)) to provide a stream of east offsets as an east elongated sequence g(i_(e)). Similarly, the north sequence elongator 1144N processes the north dither level standard deviation an with the sequence w(i_(n)) to provide a stream of north offsets as a north elongated sequence g(i_(n)).

The vertical sequence elongator 1144U processes the up dither level standard deviation au with the sequence w(i_(u)) to provide a stream of up offsets as an up elongated a sequence g(i_(u)). In an embodiment, the mean bias error vector components k_(0e), k_(0n) and k_(0u) are used by the east, north and up sequence elongators 1144E, 1144N and 1144U, respectively, to bias the mean values of the g(i_(e)), g(i_(n)) and g(i_(u)) sequences. The steam of the complex combination of east, north and up offsets is the stream of the vector offsets that is passed to the position dither processor 1106.

The position dither processor 1106 includes an east adder 1151, a north adder 1152 and an up adder 1153. The east, north and up adders 1151, 1152 and 1153 add the east, north and up offsets to east, north and up components, respectively, of the secure position 1010 to provide the east, north and up components of the rover position 1038 that is available to a user.

The equations 2-8 show the computation of the east, north and up dither levels. Once the hzPrecInput and vtprecInput (the horizontal and vertical selected precisions, respectively) are defined, it is necessary to compute the level of position dither needed. The dither variances shown in equations 2 and 3 are the differences between the actual (intrinsic) position variances achieved with a KPP position fix for the secure position 1010 and the required (selected) position variances.

hzVarDither=hzPrecInput*hzPrecInput−(Qee+Qnn)  (2)

vtVarDither=vtPrecInput*vtPrecInput−(Quu)  (3)

The Qee, Qnn, Quu are the diagonal components of the actual KPP position fix covariance matrix (typically around 1.0e-4 m² for fixed-ambiguity solutions). The hzVarDither, vtVarDither are the horizontal and vertical variances of the dither needed to provide the selected precision 1041. If the kinematic post processing is only generating a float solution, then it is possible that hzPrecOutput² is less than the current horizontal position precision (Qee+Qnn), in which cause no error is added.

The relative east/north dither level may be set based on the ratio of the east and north covariance elements in the position fix as shown in equations 4 and 5.

scaleEastDither=Qee/(Qee+Qnn)  (4)

scaleNorthDither=Qnn/(Qee+Qnn)  (5)

The east, north and up coordinates of each KPP position fix are dithered separately by east, north and up dithering processes. The dithering may use different east, north, up sequences (g_(e)(i), g_(n)(i), g_(u)(i)) or the same sequence having different epoch offsets (g(i_(e)), g(i_(n)), g(i_(u))). The standard deviations used to generate the east, north and up dither processes are then computed according to equations 6, 7 and 8.

$\begin{matrix} {\sigma_{e} = {{scaleEastDither} \times \sqrt{hzVarDither}}} & (6) \\ {\sigma_{n} = {{scaleNorthDither} \times \sqrt{hzVarDither}}} & (7) \\ {\sigma_{u} = \sqrt{vtVarDither}} & (8) \end{matrix}$

Using Satellite Geometry to Control Dither Level

An optional design of the dither process uses the prevailing satellite geometry to set the level of dithering. The solution DOP provides an unweighted measure of the satellite geometry and therefore provides an input that may be used to adapt the output statistics for the user rover position 1038. In an embodiment, the acceptable range of positional dilution of precision (PDOP) for KPP positioning is up to 7.0. In this case the selected precision 1041 would become worse as DOP increased and better as DOP decreased.

FIG. 21 is a block diagram of an optional dilution of precision (DOP) scaler 1160 that may be included in the vector offset generator 1104 to use satellite geometry to control degradation. The DOP computed by the post processor 1102 (FIG. 18) for the secure position 1010 includes a horizontal DOP (HDOP) and a vertical DOP (VDOP). The DOP scaler 1160 includes a horizontal precision selector 1162 and a vertical precision selector 1164. The horizontal precision selector 1162 determines the horizontal selected precision by multiplying the HDOP by a selected HDOP scale factor and adding a selected HDOP offset to provide the hzPrecInput (FIG. 20). Similarly, the vertical precision scaler 1164 determinants the vertical selected precision by multiplying the VDOP by a selected VDOP scale factor and adding a selected VDOP offset to provide the vtPrecInput (FIG. 20).

The selected horizontal and vertical precisions are obtained via equations 9 and 10.

hzPrecInput=HDOP_SCALE_FACTOR*HDOP+HDOP_OFFSET  (9)

vtPrecInput=VDOP_SCALE_FACTOR*VDOP+VDOP_OFFSET  (10)

where HDOP and VDOP scale factors and the HDOP and VDOP offsets define how the selected precision 1041 is adjusted according to the prevailing HDOP and VDOP.

FIG. 22 is a flow chart of a method for providing a rover position 1038 having a selected precision 1041 to a user. Steps of the method may be in the form of computer-readable instructions stored on or in or carried by a medium 1200 that may be read by a processor in a particular computer or computers for executing the steps. In a step 1202 global navigation satellite system (GNSS) observable reference data is determined by measurement and computation where the reference data includes real time kinematic (RTK) GNSS reference carrier phases. In a step 1204 GNSS observable rover data is determined by measurement and computation where the rover data includes real time kinematic (RTK) GNSS rover carrier phases.

The reference and rover data is received at a computer apparatus having kinematic post processing (KPP) capability in a step 1206. Typically, the computer apparatus is in an office in a different place than the place where the rover carrier phases were observed and the post processing is performed at a later time than the observations of the rover carrier phases. In a step 1210 the reference and rover data set is post processed with KPP to compute a secure position 1010 having an intrinsic precision 1040. The term post process refers to the data set being processed at a later time in a different machine than the reference and rover carrier phases were observed. In a step 1220 the secure position 1010 is dithered with a stream of offset vectors to compute the user-available position 1038 having the selected precision 1041 that is degraded with respect to the intrinsic precision 1040.

FIG. 23 is a flow chart of a method for generating an elongated sequence of offset vectors for dithering the secure rover position 1010. Steps of the method may be in the form of computer-readable instructions stored on or in or carried by a medium 1230 that may be read by a processor in a particular computer or computers for executing the steps. In a step 1232 the user-available precision 1041 is selected. In a step 1240 the selected precision 1041 is used to compute a dither level. In a step 1242 a sequence seed is generated for a particular data set that is to be post processed. In a step 1244 a settling time constant is selected. In a step 1246 a settling mean (bias error) is selected for the particular data set.

A Normal sequence is generated in a step 1248 with the seed. In a step 1250 the seeded Normal sequence is elongated with the selected settling time constant to generate the elongated sequence of offset vectors having the computed dither level as shown in the algorithms 1 or 1b.

FIG. 24 is a flow chart of a method for dithering a secure position 1010 to provide a user-available position 1038. Steps of the method may be in the form of computer-readable instructions stored on or in or carried by a medium 1300 that may be read by a processor in a particular computer or computers for executing the steps. In the step 1210 the secure position 1010 is computed. In a step 1302 the intrinsic precision 1040 of the secure position 1010 is determined from the covariance matrix of the position fix or some other way. In an embodiment the intrinsic precision 1040 has horizontal and vertical components and the horizontal component has east and north components. In steps 1304 and 1306 a horizontal precision and a vertical precision that are available to the user are selected.

In a step 1312 a horizontal variance is computed from the selected horizontal precision that compensates for the horizontal component of the intrinsic precision. Similarly, in a step 1314 a vertical variance is computed from the selected vertical precision that compensates for the up component of the intrinsic precision. In a step 1316 the horizontal variance is scaled according to the east intrinsic precision. Similarly, in a step 1318 the north variance is scaled according to the north intrinsic precision.

In steps 1322 and 1324 east and north standard deviations are computed from the east and north variances. In a step 1326 an up standard deviation is computed from the vertical variance. In steps 1332, 1334 and 1336 respective sequences of east, north and up offsets are computed having the east, north and up standard deviations. In an embodiment the respective Normal east, north and up sequences are generated and elongated with a selected time constant having the above described seed corresponding to a reference and rover data set. In a further embodiment in steps 1337, 1338 and 1339 the east, north and up sequences, respectively, may be biased to have non-zero average levels. In steps 1342, 1344 and 1346 the east, north and up sequences are used to dither the east, north and up components of the secure position 1010, respectively, to provide the available rover position 1038 to a user.

FIGS. 25A and 25B are flow charts of a method for selecting the user-available precision 1041 based on a dilution of precision (DOP) of the secure position 1010. Steps of the method may be in the form of computer-readable instructions stored on or in or carried by a medium 1360 that may be read by a processor in a particular computer or computers for executing the steps. In the steps 1362 and 1364 horizontal DOP (HDOP) and vertical DOP (VDOP) scale factors are selected. In steps 1366 and 1368 HDOP and VDOP offsets are selected. In steps 1372 and 1374 an HDOP and a VDOP are computed for the constellation of satellites that is used for computing the secure position 1010.

In steps 1376 and 1378 the HDOP and VDOP of the secure position satellites are scaled by the HDOP and VDOP scale factors, respectively. Similarly, in steps 1382 and 1384 the scaled HDOP and VDOP are offset by adding the HDOP and VDOP offsets, respectively.

The media discussed above, 1200, 1230, 1300 and 1360, may have a tangible form such as, but not limited to, a compact disk (CD), a digital video disk (DVD), a flash memory, a hard disk, a floppy disk, or a memory chip. The media discussed above may also be in a form of a communications medium such as, but not limited to, the internet.

Although the present invention has been described in terms of the presently preferred embodiments, it is to be understood that such disclosure is not to be interpreted as limiting. Various alterations and modifications will no doubt become apparent to those skilled in the art after having read the above disclosure. Accordingly, it is intended that the appended claims be interpreted as covering all alterations and modifications as fall within the true spirit and scope of the invention. 

1. A computer apparatus for post positioning with a selected precision, comprising: a global navigation satellite system (GNSS) post processor to post process reference GNSS carrier phases from a reference system and rover GNSS carrier phases from a rover receiver to compute a secure position for said rover receiver, said secure position not available to a user; and a random process generator to generate a sequence of offset vectors to dither said secure position according to a computed dither level to provide said selected precision for a user-available position for said rover receiver.
 2. The apparatus of claim 1, wherein: the post processor is constructed to perform said post process in a different geographical place and at a later time than said rover GNSS carrier phases are observed.
 3. The apparatus of claim 1, further comprising: a seed generator to associate a seed with a particular data set of said rover carrier phases; and wherein: the random process generator is constructed to use said seed to generate the same said sequence of said offset vectors each time said sequence of said offset vectors is generated for said particular data set.
 4. The apparatus of claim 1, wherein: the random process generator includes a Normal sequence generator to provide a Normal sequence; and an elongator to elongate said Normal sequence with a selected settling time constant.
 5. The apparatus of claim 4, wherein: the random process generator includes a bias generator to generate a bias error vector to associate with a particular data set of said rover carrier phases; and said elongator is constructed to use said bias error vector to bias a mean of said sequence of said offset vectors for said particular data set.
 6. The apparatus of claim 1, wherein: the random process generator includes east, north and up random processors constructed to use computed east, north and up components of said dither level to provide corresponding sequences of east, north and up scalar offsets, respectively, said sequence of said offset vectors comprised of said sequences of said east, north and up scalar offsets.
 7. The apparatus of claim 6, further comprising: a position dither processor including east, north and up adders to add said sequences of said east, north and up scalar offsets to east, north and up components of said secure position to provide said user-available position.
 8. The apparatus of claim 1, further comprising: a vector offset generator to compute said dither level to degrade an intrinsic precision of said secure position to provide said selected precision.
 9. The apparatus of claim 8, wherein: the vector offset generator includes an accuracy leveler to compensate said selected precision for a variation in said intrinsic precision by computing a smaller said dither level when said intrinsic precision is less precise.
 10. The apparatus of claim 9, further comprising: a DOP scaler to compute said selected precision by using a selected DOP scale factor and a selected DOP offset for processing a dilution of precision (DOP) of a constellation of satellites used for determining said secure position.
 11. A method for post positioning with a selected precision, comprising: post processing reference GNSS carrier phases from a reference system and rover GNSS carrier phases from a rover receiver for computing a secure position for said rover receiver, said secure position not available to a user; and generating a sequence of offset vectors to dither said secure position according to a computed dither level for providing said selected precision for a user-available position for said rover receiver.
 12. The method of claim 11, wherein: post processing is performed in a different geographical place and at a later time than said rover GNSS carrier phases are observed.
 13. The method of claim 11, wherein: generating said sequence includes associating a seed with a particular data set of said rover carrier phases; and using said seed for generating the same said sequence of said offset vectors each time said sequence of said offset vectors is generated for said particular data set.
 14. The method of claim 11, wherein: generating said sequence includes providing a Normal sequence; and elongating said Normal sequence with a selected settling time constant.
 15. The method of claim 14, further comprising: generating a bias error vector to associate with a particular data set of said rover carrier phases; and wherein: generating said sequence includes using said bias error vector to bias a mean of said sequence of said offset vectors for said particular data set.
 16. The method of claim 11, wherein: generating said sequence includes using computed east, north and up components of said dither level for providing sequences of east, north and up scalar offsets, respectively, said sequence of said offset vectors comprised of said sequences of said east, north and up scalar offsets.
 17. The method of claim 16, further comprising: adding said sequences of said east, north and up scalar offsets to east, north and up components of said secure position to provide said user-available position.
 18. The method of claim 11, further comprising: computing said dither level to degrade an intrinsic precision of said secure position for providing said selected precision.
 19. The method of claim 18, further comprising: compensating said selected precision for a variation in said intrinsic precision by computing a smaller said dither level when said intrinsic precision is less precise.
 20. The method in claim 19, further comprising: computing said selected precision by using a selected DOP scale factor and a selected DOP offset for processing a dilution of precision (DOP) of a constellation of satellites used for determining said secure position.
 21. A computer-readable medium having computer-executable instructions stored or carried thereby that when executed by a processor perform a method comprising steps of: post processing reference GNSS carrier phases from a reference system and rover GNSS carrier phases from a rover receiver for computing a secure position for said rover receiver, said secure position not available to a user; and generating a sequence of offset vectors to dither said secure position according to a computed dither level for providing said selected precision for a user-available position for said rover receiver.
 22. The medium of claim 21, further comprising: generating said sequence includes associating a seed with a particular data set of said rover carrier phases; and using said seed for generating the same said sequence of said offset vectors each time said sequence of said offset vectors is generated for said particular data set.
 23. The medium of claim 21, further comprising: generating said sequence includes providing a Normal sequence; and elongating said Normal sequence with a selected settling time constant.
 24. The medium of claim 23, further comprising: generating a bias error vector to associate with a particular data set of said rover carrier phases; and wherein: generating said sequence includes using said bias error vector to bias a mean of said sequence of said offset vectors for said particular data set.
 25. The medium of claim 21, wherein: generating said sequence includes using computed east, north and up components of said dither level for providing sequences of east, north and up scalar offsets, respectively, said sequence of said offset vectors comprised of said sequences of said east, north and up scalar offsets. 